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The metric structure of bosonic scale-free networks and fermionic Cayley-tree networks is an- 
alyzed focousing on the directed distance of nodes from the origin. The topology of the netwoks 
strongly depends on the dynamical parameter T, called temperature. At T = oo we show ana- 
lytically that the two networks have a similar behavior: the distance of a generic node from the 
origin of the network scales as the logarithm of the number of nodes in the network. At T = 
the two networks have an opposite behavior: the bosonic network remains very clusterized (the 
distance from the origin remains constant as the network increases the number of nodes) while the 
fermionic network grows following a single branch of the tree and the distance from the origin goes 
as a power-law of the number of nodes in the network. 



PACS numbers: 89.75.-k, 89.75.Hc 

Complex networks representing systems of interact- 
ing units have been studied and classified [1,2] according 
to their different geometrical and topological properties. 
In between complex networks scale-free networks, with 
power-law connectivity distribution, have been found to 
describe different systems of nature and society [3-5] . Re- 
cently several models [3-6] have been formulated that 
generate such structures, the prototype of them being 
the B A model [7] . Scale-free networks are particulary in- 
teresting because their highly inhomogeneous structure 
induces peculiar effects in the dynamical models that 
can be defined on them, as the absence of percolation 
[8] and epidemic threshold [9] ,the infinite Curie tempera- 
ture for the paramagnetic to ferromagnetic transition in 
the frame of the Ising model [10-13] and the good asso- 
ciative memory of the Hopfiled model defined on network 
with large average connectivity [14]. Besides, the investi- 
gation of the metric stucture [15-19] of these networks is 
of great interest. It was first empirically found [15] and 
then analytically derived [16,17] that scale-free networks 
are characterized by having a mean distance < d > be- 
tween nodes scaling like the logarithm of the system size 
iV, < d >- \og{N). 

The similarities between the structure of a BA net- 
work and a traditional Cayley tree have been recently 
studied. It has been found that they share high similar- 
ities. In fact they are both generated by the subsequent 
addition of a same elementary unit (a node connected to 
m links) attached to the rest of the network in the direct 
or in the reversed direction [20] . The symmetry between 
these two types of networks is evident if we assume that 
each node i has an innate quality or "energy" and 
that the dynamics is parametrized by a variable T called 
temperature that introduce a "thermal noise" as it has 
been done in self-organized models [21-23]. It is possible, 
then to observe that the BA model becomes a limiting 
case of a scale-free network described by a Bose distri- 
bution of the energies to which the incoming links point 
[24] , while the growing Cayley tree network is described 
by a Fermi distribution of the energies at the interface 
[25]. Quantum networks (the bosonic scale-free network 
and the fermionic Cayley tree network) evolve around a 



well defined core of initial nodes. 

In this paper we focus our attention on the metric 
structure of quantum networks and in particular on their 
size, i.e. we extimated the distance measured over di- 
rected paths of a generic node i from the origin of the 
network and its dependence on the time ti in which the 
node i has been added to the network. We derive the 
expression for the average value of the distance < £{ti) > 
of node i introduced at time ti from the origin, mesuread 
on directed paths, in quantum networks at T = oo. We 
find in agreement with [15-17] that < £{ti) >~ \og{ti) in 
the bosonic scale-free network and we show for T = oo a 
similar behavior in fermionic networks. 

At different values of T the topology of the network 
change drastically. For energy distribution functions 
p{e) — > for e — > the exponent 7 of the power-law 
connectivity distribution P{k) ~ goes from 7 = 3 at 
T = 00 to 7 = 2 for T = and there is a phase transition 
at a cirtical temperature Tc below which a finite fraction 
of all the links is connected to a single node. Assuming 
that the bosonic network is a reasonable model for grow- 
ing scale-free networks, we will have, for example, that 
the citation network [26] with 7^3 would correspond to 
aT = 00 dynamics while the incoming component of the 
World- Wide- Web with 7^^ = 2.1 [7,27] would correspond 
to a low temperature dynamics. As T decreses also the 
behavior of < £{ti) > changes and the distance of node 
i from the origin depends less strongly on the time ti 
of its arrival. At sufficiently low temperature i{ti) re- 
mains constant as a function of ti. On the contrary in 
the fermionic network, at T = 0, when the dynamics be- 
comes extremal, the network evolves far away from the 
origin and the distance of a node i from the origin of the 
network grows as a power-law of the time ti of its arrival 
in the network. 



A. Distances from the origin in the bosonic network 

The bosonic network is a generalization of the well 
known BA network [7]. In this model, a new node with 
m links is added to the network at each time step. Each 
node i has an innate quality or 'energy' Ci extracted 
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from a probability distribution p(e). The way the new 
links arc attached follows a generalized preferential at- 
tachment rule: the probability Hj that an existing node 
j acquires a new link depends both on its connectivity 
kj{t) and its energy ej, i.e. 



(1) 



with the parameter (3—1 /T tuning the relevance of the 
energy €j with respect to the connectivity kj{t). This 
network displays a power law connectivity distribution 
P{k) ~ k~'^ with 7 G [2, 3] depending on the p(e) distri- 
bution and the inverse temperature (3 [24] . In the T = oo 
limit (/3 = 0) the network reduces to a scale-free BA net- 
work [7] with an average connectivity fcj of node i that 
grows in time as a power-law with exponent 1 /2 



ki{t) 



(2) 



where U is the time in which node i was added to the 
network. The probability pij that two nodes i and j are 
connected by a hnk can be calculated from Eq. (1). At 
T = 00 (/3 = 0), taking into account that at each time m 
new links are added to the network, the probability Pij 
is given by m times Eq. (1). After substituting (2) into 
Eq. (1) we obtain, 



This means that the mean distance between node i and 
a node io, such that =1, calculated only on directed 
paths, follows a Poisson distribution with average size 
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FIG. 1. Mean distance from the origin < i(ti) > of the 
nodes arrived at time U in a BA network (a bosonic network 
at T = oo) with m — 1,2. The solid lines indicate the theo- 
retical prediction Eq. (6). 
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The number of directed paths of length £ connect- 
ing node i, introduced in the network at time ti, to 
node io belonging to the original core of the network 
{tio = 1) is given by the mean value of the number of 
paths connecting node i to io and passing through the 
points io, ii,i2f ■ ■ = i with tn+i > tn- Indicating 

each node with the time of its arrival in the network, the 
probability of any directed path is given by the product 
n^=iPri-i,n with tn-i < tn- In Order to find n£{ti) we 
should sum over all possible paths. Replacing the sum 
over the nodes with the integrals over tn, we obtain for 
ne{ti), 



ne{ti) = dti dt2.-. dte-iPQ,iPi 
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Using ^0 = 1 and Eq. (3) valid in the T = oo limit, we 
obtain 
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The distribution of the number of directed paths of 
length £ starting from the origin of the network is pro- 
portional to the integral of ne{ti) over ti 



Pi£)o. rn,it')dt' = {-mY(l-^-^g^]. (7) 



r(^) ; 



In Fig. 2 we show the agreement between the numerical 
results and Eq. (7) for a bosonic network of 10^ nodes at 
T = oo and m = 1, 2. 




FIG. 2. Distribution of the number of directed paths of 
length i in a BA network (a bosonic network at T = oo) for 
m = 1,2. The solid lines are the analytic predictions Eq. (7). 
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The topology of the bosonic network changes as a func- 
tion of the temperature T = 1 / (3. In particular for a dis- 
tribution p{e) such that p(e) — > as e ^ we know [24] 
that there is a critical temperature below which the 
network has a topological transition and its structure is 
dominated by a single node that grabs a finite fraction 
of all the hnks. For T > Tc {(3 < 13c) the network is in 
the so-called "fit-gets-rich" phase (FGR phase) while for 
T < Tc {(3 > j3c) the network is in the so-called Bosc- 
Einstein condensate phase (BE phase). To visualize this 
transition we have plotted the bosonic network with an 
energy distribution 

P^') = ^/^ ^^(O'l) (8) 

where 6 = 0.5, m = 1, above (Fig. 3a) and below (Fig. 
3b) the phase transition. 



(a) 



(b) 




FIG. 3. Graphic representation of the bosonic network. 
Graph (a) represents a network with m = 1, energy distri- 
bution given by Eq. (8) where 6 = 0.5 at T = oo {P = 
-the BA network) and graph (b) represents a network with 
the same parameters m and 9 but with T = 1//3 = 0.33 (net- 
work in the Bosc-Einstein condensate phase). The number of 
nodes both networks is A'^ = 10^. 

The network has been designed in order to underline 
the hierarchical structure of the network. Starting from 
the single node at the origin of the tree, we have placed 



all the nodes that are directly attached to it on a semicir- 
cle of unitary radius, each node i separated from the next 
one by the angle Aa^ proportional to its connectivity, i.e. 

Aa. = (9) 

where N{i) are the nearest neighbors of node i added at 

a time tj > ti. We have repeated the same construction 
for all the nodes of the network in such a way that all the 
nearest neighbors of node i are on a semicircle of radius 
with 

n = rk ^ . (10) 

where k is the node to which the node i has been at- 
tached at time ti. Prom Fig. 3 it is clear the change in 
the topology of the network at the critical temperature 
with the emergence of a single node that grabs a finite 
fraction of all the links in the Bose-Einstein condensate 
phase. 

As the topology of the network changes, the behavior 
of as a function of ti changes too. In fact we have 

< £(i) >= a(/3) log(ii) (11) 

where the coefficient a{(3) is a decreasing function of the 
inverse temperature (3 = 1/T. 
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FIG. 4. Distances form the origin in a bosonic network with 
m = 2 and 9 = 0.5 as a function of (3. 

In Fig. 4 we report < £{i) > for a bosonic network with 
p{e) of the type (8) with = 0.5 and m = 2 at differ- 
ent values of the inverse temperature f3, above and below 
the critical value 0c = 1-7 [24]. In oredr to illustrate the 
changemont of the topology of the network above and 
below the critical temperature in Fig. 5 we report the 
behavior of different relevant structural quantities for a 
network of size N = IC*. 9 = 0.5 and with m = 2 around 
the critical inverse temperature temperature /?c = 1-7 . 
We report the fraction of links attached to the most con- 
nected node, the exponent of the power-law component 
of the connectivity distribution, the clustering coefficient 
and the coefficient a(/3). The fraction of links attached 
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to the most connected network kmax{P)/N is the order 
parameter of the FGR-BE phase transition and increases 
as a function of (3. The data reported in Fig. 5 are medi- 
ated over 100 runs. The connectivity distribution of the 
bosonic network contains a power-law component plus 
a pick indicating the condensation phenomena that ap- 
pears for /3 > Pc- In Fig- 5 we report the exponent 7(/3) 
of the power-law component of the connectivity distribu- 
tion that decreases as a function of /? with an asymptotic 
value of 7 = 2. The data reported on Fig. 5 are mediated 
over 100 runs. The clustering coefficient C{0) increases 
at the transition point while the coefficient a{/3) slowly 
decreases saturating toward a zero value for (3 > 3.0. The 
data reported in Fig. 5 for these last two quantities data 
are mediated over 10 run. 
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where the sum in the denominator is extended to all 
nodes j at the interface Int{t) at time t. The model de- 
pends on the inverse temperature /? = 1/T. In the /3 ^ 
limit, high and low energy nodes are equally probable to 
grow and the model reduces to the Eden model. In the 
(3 ^ oo limit the dynamics becomes extremal and only 
the nodes with the highest energy value are allowed to 
grow. In this case the model reduces to invasion perco- 
lation [28,29] on a Cayley tree. 

Let us assume that /3 = (T = oo). The probability 
Hi that a node i of the interface Int{t) grows at time t 
is given by 



n, 



1 



Nintity 



(13) 



where Nint{t) is the total number of active nodes. Since 
at each timestep a node of the interface branches, and 
m new nodes are generated, after t timesteps the model 
generates an interface of Njnt{t) nodes, with 



Njutit) = (m - l)t + 1. 



(14) 



Lets denote by p{t, ti) the probability that a node, born 
at time U is still at the interface at time t. Since every 

node i of the network grows with probability 11^ Eq. (13) 
if it is at the interface, in mean field p(t,ti) follows 



dp{t, tj 
dt 
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Replacing (14) in (15) in the limit t 
solution 



(15) 



oo we get the 



FIG. 5. Relevant structural quantities in a bosonic network 
of size N = l(f with m = 2 and d = 0.5 as a function of the 
inverse temperature (3. 



B. Distances from the origin in the fermionic 
network 

The fermionic network [25] is a growing Cayley tree, 
where the innate qualities of the nodes (energies) define 
their different branching tendency. Starting at time t = 1 
from a node iq at the origin of the network, the node iq 
at time t = 2 grows and m new nodes are directly con- 
nected to it. Each node i has an energy extracted 
from a given p(e) distribution. At each timestep a new 
node with connectivity one (at the interface) is chosen to 
branch, giving rise to m, new nodes. 

We assume that nodes with higher energy are more 
likely to grow than lower energy ones. In particular we 
take Ilj, the probability that a node i of the interface 
(with energy e^) grows at time t, to be 



^jeint(t) 



(12) 



p{t,ti) = 



l/(m-l) 



(16) 



Consequently each node i that arrives at the interface at 
time ti, remains at the interface with a probability that 
decreases in time as a power-law. The probability pij 
that a node i is attached to a node j (arrived in the net- 
work at a later time tj > ti) is given by the r. h. s. of Eq. 
(13) calculated at time tj. Taking into account Eq.(16) 
and the rate m of the addition of new nodes, we obtain 
for pij 



Pi,j = 



m 



(to - l)tj 



l/(m-l) 



(17) 



The number n^(ti) of paths of length i that con- 
nect a node i, introduced at time ti, to the origin io 
is given by the average number of paths connecting a 
node i to a node io and passing through the points 
io,i\i2, ■ ■ - ii-i^it = i with tn+i > tn- Indicating each 
node with the time of its arrival in the network and the 
sum over the nodes with the integrals over tn, we obtain 
for ne{ti). 
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ni{ti)= I dti dt2... dte_iPoiPi,2---Pe-i,e 

Jl Jti Jte-2 

(18) 



and, using (17) we obtain, with a calculation analog to 

(5), 




t. 



FIG. 6. Distance from the origin in a fermionic network at 
T = oo {13 — 0) for networks oi N = 10* nodes with p{e) uni- 
form between zero and one and with m — 2,3, 4. The soUd 
Unes are the theoretical predictions Eq.(20). 

This means that the mean distance between a node U 
and the origin follows a Poisson distribution with average 
size 



<i{ti)>=^^\og{ti). (20) 



As in the bosonic network, in the fermionic network at in- 
finite temperature (/3 = 0) the distance < i{ti) > of node 
i from the origin grows logarithmically with the time U. 

In Fig. 6 we report the analytical simulation of a 
fermionic network with p{e) uniform between zero and 
one, at T = oo and m = 2, 3, 4. 

As the temperature decreases, the topology of the net- 
work changes drastically, in Fig. 7 wc show the Caylcy 
tree with p{e) = 1, e € (0, 1) and to = 2 at infinite tem- 
perature (/3 = 0.0) and at low temperature (/3 = 20.0). 
At high temperature the network grows homogeneously 
in each direction while at low temperature it evolves fol- 
lowing only a single branch of the tree. 




FIG. 7. The fermionic network with m = 2 at T = oo 
(Graph (a)) and at temperature T = 0.05 (Graph (b)).The 
number of nodes in both networks are N = 10*. 

The distance of a node i from the origin of the net- 
works grows logarithmically with at T = oo (/3 = 0). 
As the temperature decreases the behavior of < (.{ti) > 
gets steeper. For p(e) uniformly distributed between zero 
and one, in the extremal case T = (/3 = oo) when the 
node of highest energy grows deterministically at each 
time step, we have a dramatic change in the behavior 
and < £(i) > grows as a power-law of ti, 

< e{i) >oc {ti)i (21) 

C = 0.55 ± 0.05 from the numerical results reported in 
Fig. 8. 

In conclusion we have shown that bosonic and 
fermionic network are not only simmetricaly built [20] 
but also at T = 00 they are characterized by a distance 
< (.{ti) > from the origin that grows like the logarithm of 
the time tt. On the contrary, in the limit T — they be- 
have in a opposite way: the bosonic network stays highly 
clusterized with a distance from the origin that remains 
constant as the network evolves, in the fermionic network 
the distance < ((ti) > grows like a power-law of the time 
ti. 
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FIG. 8. Distance from the origin in a femionic network with 
m = 2 and p(e) uniformly distributed between zero and one 
at different temperatures. At l3 — {T — oo) we have the 
predicted logarithmic behavior Eq. (20) while in the extremal 
case P — oo {T = 0) the network grows as a power-law of the 
network size. The solid line in the Inset is the power-law fit 
Eq. (21) with C = 0.55 ± 0.05. 
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